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ABSTRACT 

Context. Recent 1.2-mm continuum observations have shown the giant Hn region NGC3576 to be embedded in the centre of an 
extended filamentary dust-cloud. The bulk of the filament away from the Hn region contains a number of clumps seen only at 
(sub-)millimetre wavelengths. Infrared and radio observations of the central star cluster have uncovered evidence of sequential star- 
formation leading us to believe that the adjacent clumps may host massive protostellar objects at a very early stage of evolution. 
Aims. We have investigated the physical and chemical conditions in the dusty clumps with the goal of characterising their star-forming 
content. 

Methods. We have used the Australia Telescope Compact Array (ATCA) to image the cloud for the NH3 (1,1), (2,2) and (4,4) tran- 
sitions, 22 GHz H 2 masers, and 23 GHz continuum emission. The 70-m Tidbinbilla dish was used to estimate the total integrated 
intensity of NH 3 . We also utilised the 22-m Mopra antenna to map the region for the molecular lines 13 CO(1-0), C 18 O(l-0), 
HCO + (1-0), H 13 CO + (1-0), CS (1-0) and N 2 H + (1 -0). 

Results. Emission from dense molecular gas follows the morphology of the 1.2-mm dust emission, except towards the central ionised 
region. The Hn region is observed to be expanding into the molecular cloud, sweeping up a clumpy shell of gas, while the central 
star cluster is dispersing the molecular gas to the east. Analysis of the NH 3 data indicates that temperature and linewidth gradients 
exist in the western arm of the filament. Temperatures are highest adjacent to the central Hn region, indicating that the embedded 
cluster of young stars there is heating the gas. Six new H2O masers were detected in the arms of the filament, all associated with 
NH3 emission peaks, confirming that star-formation has begun within these cores. Core masses range from 5 to 516M Q and most 
appear to be gravitationally bound. Complementary results by Andre et al. (2008) imply that seven cores will go on to form massive 
stars between 15 and 50 M Q . The large scale velocity structure of the filament is smooth, but at least one clump shows the signature of 
inward gas motions via asymmetries in the NH3 (1,1) line profiles. The same clump exhibits an enhanced abundance of N 2 H + , which 
coupled with an absence of CO indicates depletion onto the dust grain surface. 

Conclusions. The Hn region at the heart of NGC 3576 is potentially triggering the formation of massive stars in the bulk of the 
associated cloud. 
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1. Introduction 

1.1. Triggered star-formation 

A necessary precursor to massive star for mation is the exis tence 
of dense clumps of self-gravitating gas. Elmegreen (1998]) de- 
veloped the first coherent picture of sequential star formation, in 
which the formation and collapse of these clumps within giant 
molecular clouds is triggered by an external event. Three distinct 
triggering mechanism are considered: 

1 . Globule squeezing: compression of pre-existing clumps, e.g. 
due to a propagating shock-wave from a supernova. 

2. Cloud-cloud collisions: two molecular clouds collide result- 
ing in gravitational instabilities. 

3. Collect and collapse: accumulation of gas into a shell or 
ridge, and subsequent fragmentation and collapse. 

Examples of the third mechanism are observed to occur on 
the edges of Hn regions. Created by the far-UV radiation from 
a young OB-cluster, the hot ionised gas expands rapidly, sweep- 
ing up a shell of dense molecular material before it. A photon 
dominated region (PDR) is created at the interface between the 
Hii re gion and the molecular cloud (see Hollenbach & T ielensl 
1999 for a review of PDRs). If the PDR is over-pressured com- 
pared to the bulk of the molecular gas, shocks are driven into 
the dense neutral medium of the c loud, potentially lea ding to its 
fragmentation into clumps (e.g. Urqu hart et al.1 12006 and refer- 
ences therein). Subsequent disturbances may then trigger col- 
lapse and hence star formation. Material swept into a clump 
shields the column of molecular material behind it, forming 
giant pillars, the best known ex ample of which are the 'ele- 
phant trunks' in the Eagle Nebula (I White et alj 19991: [Allen et alj 
1 19991) . Some observed exa mples of swept-up she lls around 
Hii regions include Sh2 -212 dDeharveng et alj|2008l) . RCW120 
dDeharveng etaLjj2009b . W5 dKarr & MartirJ 12001. R CW49 
dWhitnev et al.ll2004l) . and RCW79 dZavagno et al.ll2006l) . 

It is important to note that energy injected into molecular 
clouds by newly formed stars may instead quench the star for- 
mation process. Supersonic particle-winds from massive stars 
account for ~ 0.1 - 1.0 per cent of the stellar luminosity an d act 
to disperse molecular material and dust (e. g„ iGenzelll 199 ID . The 
two most important feedback mechanisms, photo-ionisation and 
cluster winds, compete with each other to disrupt the cloud. On 
the scale of giant molecular clouds, feedback from newly formed 
stars is responsible for regulating the star-formation rate and 
hence the evolution of Galactic structure. 

In this paper we present new multi-wavelength observations 
of the giant Hii region NGC 3576, which is thought to be under- 
going sequential star formation (see Section [TTZl l. We examine 
the evidence for star formation in the nearby molecular environ- 
ment. 

1.2. Prior observations of NGC 3576 

The environment of the giant Hii region NGC 3576 (also known 
as RCW57, G291.3-0.7 and IRAS 1 1097-6102) has been ex- 
tensively studied at infrared and radio wavelengths. Figure [Tj-a 
presents the 5.8//m SPITZER IRAC image of the whole star- 
forming complex with the Hii region at the centre. The hour- 
glass shaped emission feature running north-south across the 
centre of the image likely correspond to a bipolar cavity being 
evacuated by the central embedded cluster. The giant molecular 
cloud which hosts NGC3576 manifests itself as large extincted 
areas to the north-east and south-west, and is at a distance of 



2.4 kpc dPersi et alj|1994l) . The extent of the Hii region is il- 
lustrated by the thick black contour, corresponding to the 3-cr 
(5 4mJy/beam) level in the 3.4-cm free-free emission mapped 
bv lde Pree et al.1 d 19991) . FigureQJb shows the distribution of the 
3.4-cm emission plotted over the 2.17/mi infrared image from 
2MASS. The ionised gas covers the brightest infrared emission, 
peaking sh arply in the west a n d exte nding ~ 2' to the north-east 
and south. iFrogel & Perssonl d 19741) discovered five near- and 
mid-infrared sources towards the ionised gas (plotted as circles 
in Figure Q]), with the brightest source (IRS 1) located adjacent 
to the main peak of the radio emission . Furth er high resolution 
photometric observations bv lPersi et alj d 19941) have revealed the 
presence of a very young, deeply embedded cluster with 130 
members in the same region. Spectra o f these sources exhibit 
deep 9.7/im silicate absorption features dMoorwood & Salinaril 
1981), leading to their interpretation as pre-main-sequence ob- 
jects. 

There is contr oversial evidence for sequential star forma- 
tion in NGC 3576. |Persi et alj d 19941) showed that a steep near- 
IR colour gradient exists in the embedded cluster, with the 
reddest, most deeply embedded sources in the west. This im- 
plies that star-formation began in the east, gradually moving to 
the location of the strong radio peak in the west. Such an in- 
terpretation is supported by the exist ence of an electron tem- 
perature gradient in the ionised gas dde Pree et alj [l999) indi- 
cating that the youngest and hottest stars are locate d near the 
sharp western edge of the Hii region dHiellming|1966l) . High he- 
lium abundance in the east may also indicate the presence of a 
popul ation of older stars and their mass loss ( Hanson & Contil 
U993I) . Bright 22 GHz water masers, commonly foun d in out- 
flows , have b een found adj a cent to the main radio peak dCaswelll 
120041) . while iNorris et alj d 19931) also detected two 6.67 GHz 
methanol maser sites near IRS 1 and IRS 3. Shocked molecu- 
lar hydrogen emission was measured peaking towards IRS 1 by 
Oliva & Moorwoodld 19861) . also pointing to the existenc e of out- 
flows. IRS 1 was later resolved in to three components bv lMonetH 
d 19921) and lBarbosaetaT] d2003l) . one of which is seen through 
scattered light from a cavity or a dusty disk dMonetil 19921) . None 
of the detected infrared sour ces contributes sign ificantly to the 
ionisation of the Hii region, dBarbosa et al] |2003l). however, re- 
cent Chandra observations dTowns lev 2006) have revealed mul- 
tiple deeply embedded hard X-ray sources, which may provide 
the extra ionisation needed. 

The Hii region was mapped at 1 -mm w avelengths for contin- 
uum emission, with a resolution of 1', bv lCheung et alj (jj980), 
who found that the flux density ratio compared to 40-350/im was 
consis tent with optical ly thin thermal dust emission. More re- 
centlv. lHill et alJd2005l) used the SIMBA0 bolometer to map 1.2- 
mm continuum emission at a resolution of ~ 24", as shown in 
Figure [TJ-c and via white contours on Figure [TJ-a. The SIMBA 
field covers an area of 4'x6', centred on the radio peak and 
reveals that the Hii region is embedded in a filamentary struc- 
ture of cool dust, running north-east to south-west. The Hii re- 
gion is prominent at the centre, however, several bright knots are 
apparent along the length of the filament. These clumps, des- 
ignated S 1 - S5 in this work, fall along the dark lane travers- 
ing the nebulosity in the near-infrared image, and are coincident 
with infrared-dark clouds (IRDCs) in the mid-infrared (e.g., the 
5.8ytzm IRAC image). IRDCs have been found throughout the 
Galactic plane with masses in excess of 30 M and it has been 
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Fig. 1. Overview of the NGC 3576 star forming region, (a) SPITZER 5.8pm image showing the extent of the 3.4-cm radio emission 
(black contour, ~3 percent) and the 1.2-mm continuum emission (white contour, ~5 percent). CH3OH and H2O maser positions are 
m arked with 'x' and ' +' sym bols, respectively. Circles mark the positions assigned to the original infrared sources (IRS) discovered 
by iFrogel & Perssonl d 19741) . while a ' □ ' marks the peak of the radio-contin uum emission, (b) C ontours of the 3.4-cm free-free 
emission from ionised gas in the giant Hn region at the heart of th e complex dde Pree et alj | 1999h plotted over the 2.2jum K-band 
image from 2MASS. (c) 1.2-mm continuum image of the complex dHill et al.ll2005l) showing the dusty filament running across the 
Hn region. Dashed circles indicate the five fields observed with the ATCA, centred on the 1.2-mm clumps. 



sugg ested they are the cold precursors to m assive star clusters 
(e.g JPillai et1d1l2006llR7thborne et al.ll2006l) . 

NGC 3576 constitutes an ideal laboratory in which to study 
the process of massive star formation. In particular, we aim to 
investigate if the Hn region is interacting with the dusty filament 
and if star-formation has been triggered in the dense knots along 
its length. Here we present new observations of the whole fila- 
ment in several molecular tracers, with the goal of constraining 
the star-forming properties of these new clumps (SI, S2, S4 & 
S5), and determining the effect of the Hn region on the complex. 



2. Observations and data reduction 

Data from three telescopes, the Australia Telescope Compact 
Array (ATCA), the 22-m Mopra telescope and the 70-m 
Tidbinbilla telescope, were combined to assemble a picture 
of NGC 3576. The ATCA was used to obtain high resolution 



NH3 maps as a probe of the density and temperature structure 
of the clumps within the filament. We simultaneously searched 
for UCHn regions via their 23 GHz free-free continuum emis- 
sion. We used the Tidbinbilla telescope to map the extended 
NH3 emission, providing a measure of the total column of 
NH3 in the filament and an estimation of the 'missing flux' in 
the interferometer images. We utilised the Mopra telescope to 
map the complex in the 3-mm lines 13 CO(1-0), C 18 O(l-0), 
HCO + (1-0), H 13 CO + (1-0), N 2 H + (l-0) and CS (2-1), with 
the goal of probing the physical and chemical conditions in each 
of the dusty clumps. Line rest frequencies and electronic con- 
stants for each transition are noted in Table lA.ll in the Appendix. 

2.1. ATCA observations 

Observations were made with the ATCA in three blocks dur- 
ing the years 2003-2005. Table Q] summarises details of the 
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Table 1. Details of the Australia Telescope Compact Array observing runs. 



UT 


Array 


Beam" 


Transitions 


Frequency 


Bandwidth 


No. 


K/(Jy/bm) 


Date 


Config. 


Size (") 




(GHz) 


(MHz) 


Channels 




2003 Aug 28 


EW367 


7.6x6.2 


NH 3 (U) 


23.694480 


8 


512 


46.02" 








Continuum 


23.694500 


128 


32 










H 2 (6-5) 


22.235120 


16 


16 




2004 Jul 02 


750D 


6.6x2.4 


NH 3 (4,4) 


24.139417 


8 


512 


134.47 








Continuum 


22.000000 


128 


32 




2005 Jul 17 


H75 


26.8 x 20.9 


NH 3 (1,1) 


23.694480 


8 


512 


3.89 








NH 3 (2,2) 


23.722630 


8 


128 


3.88 



" Beam size excluding the 3 km+ baselines to Antenna 6. 

13 Jy to K conversion factor for the combined EW367 and H75 data is 17.58 K/(Jy/bm). 



Table 2. Centre coordinates of the ATCA fields corresponding 
approximately to the peak 1.2-mm emission. 



Field 


RA 


DEC 




(J2000) 


(J2000) 


SI 


11" 12'" 17.00' 


-61" 17'" 40.0 s 


S2 


11"12'"04.90 S 


-61 rf 18 m 20.0 s 


S3 


ll"ll m 5 1.50 s 


-61 rf 18" , 52.0 s 


S4 


ll' , ll m 38.33 s 


-61 rf 19 m 53.9 s 


S5 


ll*ll m 33.86 s 


-61 £ '21 m 20.9 s 



dates, array-configurations and frequencies used. NH3 (1,1) data 
was obtained using the EW367 and H75 array configurations in 
August 2003 and July 2005, respectively. Two orthogonal lin- 
ear polarisations were observed at each frequency. The raw data 
were processed in a correlator, which was configured to deliver 
a bandwidth of ~ lOlkms -1 split into ~0.2kms _1 wide chan- 
nels. NH3 (2,2) was observed simultaneously with NH3 (1,1) on 
the H75 array in July 2005. Due to technical constraints only 
128 channels were available on the correlator, yielding a spec- 
tral resolution of ~ 0.8 km s over the 101 km s _1 bandwidth. In 
August 2003 the region was also mapped for 23 GHz continuum 
and 22 GHz water masers using the EW367 array. In July 2004, 
the 750D array was used to simultaneously map the NH3 (4,4) 
transition and 23 GHz continuum. The correlator configuration 
was identical to that used to obtain the NH3 (1,1) data. 

The primary beam of the ATCA at 23 GHz has a full width 
half maximum (FWHM) of 2.5'. In order to cover the extent 
of the 1.2mm continuum emission we observed five overlap- 
ping fields, marked by dashed circles in Figure [T}d and centred 
approximately on the coordinates of the dusty SIMBA clumps. 
Table|2]lists the coordinates of the pointing centres. Each of the 
five fields was observed for 10 minutes in turn, over the course of 
one hour and this pattern was repeated for 10 hours, giving a total 
of 1 .7 hours on each position. In order to correct for fluctuations 
in the phase and amplitude caused by atmospheric and instru- 
mental effects, a strong phase calibrator was observed for two 
minutes before and after changing fields. The instrumental con- 
tribution to the bandpass shape was measured by integrating on a 
strong continuum source (e.g., 0420-014) and was subsequently 
subtracted from all spectra. A primary flux calibrator (1934-638 
or Uranus) was observed once per observation period, to allow 
the absolute calibration of the flux scale. 

T he data were reduced using the miriad package (ISault et al.l 
Il995l) following standard procedures. During the data reduction 
the sources were assumed to be unpolarised and both polari- 
sations were averaged together. Continuum emission was sub- 
tracted from the spectral line data by using the task uvlin to fit a 



polynomial to the line free-channels. An image of the continuum 
emission was produced by combining the line-free channels in 
the NH3 data with the dedicated wide-band 22 and 23 GHz con- 
tinuum observations. Images were made using the invert task and 
natural weighting was used to minimise the noise in the image- 
plane. At this stage the ~ 3 km baselines to antenna CA06 were 
discarded as the phases were found to be decorrelated, degrad- 
ing the image quality. All images were deconvolved using the 
standard clean algorithm and, if a sufficiently bright source was 
present, several iterations of the selfcal task were applied. The 
above procedure was performed on all fields, before using the 
task linmos to mosaic the data into a single map. Finally the 
NH3 data were converted to a brightness temperature scale in 
Kelvin, using the Jy/K scaling factors noted in Table Q] 



2.2. Mopra observations 

The Mopra observations were conducted in 'on the fly' (OTF) 
mapping mode, between the months June - September, during 
2004 and 2005. An image was built up by combining overlap- 
ping scan rows, each containing 30-40 spectra. Scan rows were 
offset by half of the beam-FWHM and the scanning speed was 
slow enough to ensure Nyquist sampling in the scan direction. 

NGC 3576 was divided into an overlapping mosaic of three 
or more 5' x5' fields, positioned to cover the 1.2-mm contin- 
uum emission. The pointing centres of individual fields varied, 
depending on where molecular emission was detected, and the 
final maps were assembled by co-adding individual fields to- 
gether. Each field took ~ 80 minutes to complete, plus a further 
~ 10 minutes for pointing checks and calibration measurements. 

The signal from the receiver was processed in a digital auto- 
correlator, configured to have a bandwidth of 64 MHz divided 
into 1024 channels, which provided a velocity resolution of 
~0.2kms~' over a usable bandwidth of ~ 120kms _1 . The cen- 
tral frequency was chosen so that channel 512 was centred on 
the systemic velocity of NGC 3576 at -24 km s _1 . Observations 
were made in dual orthogonal linear polarisation mode and the 
polarisations were averaged together during the reduction pro- 
cedure. The pointing accuracy was checked using a nearby SiO 
maser before observing each field and was estimated to be better 
than 8". Calibration to the T* A scale was achieved by measuring 
the emission from a singl e hot load placed i n fro nt of the receiver 
every 20 minutes (see lLadd et all 120051 and iKutner & Ulichl 
1981). The maps were further calibrated onto the main-beam 
brightness temperature scale (Tmb) by dividing the T^ pixel val- 
ues by the mai n beam efficiency r/ m \, at the observing frequency 
(see Table 4 of lLadd et al.ll2005]) . 
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The data were reduced using the livedata and gridzilla 
packages, available from the ATNrfl livedata performed band- 
pass calibration by subtracting the preceding OFF spectrum from 
the SCAN spectra in each row. A low-order polynomial was then 
fit to line-free channels and subtracted, resulting in a smooth 
baseline at zero Kelvin. The spectra were assigned individual 
position stamps and regridded to the LSR-K reference frame, be- 
fore being written to disk. The gridzilla package was then used 
to resample the maps to a regular pixel scale, weighted accord- 
ing to the system temperature (T sys ). To grid the data we used 
a pixel-size of 6x6 arcseconds and a Gaussian smoothing ker- 
nel with a FWHM of 18", truncated at an angular offset of 36". 
The final data cubes were smoothed to an angular resolution of 
-40". 

2.3. Tidbinbilla observations 

The 70-m antenna located at the NASA Tidbinbilla Deep Space 
Communications Complex near Canberra has a limited time de- 
voted to radio astronomy. When Mars is above the horizon the 
70-m telescope is almost entirely utilised supporting the various 
NASA missions to that planet and as such has been little used by 
the radio astronomy community. In 2005 we performed the first 
frequency-switched spectral line mapping observations using the 
23 GHz K-band receiver. Fast frequency switching also provides 
better cancellation of atmospheric emission fluctuations, how- 
ever, the response of the receiver may vary considerably with 
frequency, resulting in poor spectroscopic baselines. 

We targeted the same five fields in NGC 3576 as were ob- 
served with the ATCA. The OTF method had not yet been im- 
plemented on the 70-m antenna so we observed each field as 
a square grid of 5x5 positions, spaced by half the 45" beam- 
FWHM. Each position was observed for 10 minutes and the 
pointing accuracy was checked by observing an unresolved 
planet (Jupiter) every hour. Pointing errors were usually below 
10" or 1/4 beam. Individual maps took 5-6 hours to complete, 
including 15 minutes to measure the atmospheric opacity using 
the 'skydip' method. During the observations the data were cal- 
ibrated against an ambient load (a noise diode) and the zenith 
system temperature at 23 GHz was typically 40-50K during 
the observations. We estimate the T^ flux scale is uncertain by 
< 10 per cent. The peak aperture efficiency wa s measured as 
0.48 + 5% at 22.2 GHz by iGreenhill etall {2003) and we adopt 
that value here to correct the data onto the main beam brightness 
temperature flux scale (Tmb)- To convert to units of Jy/beam the 
data were multiplied by 0.95 1 (Jy/beam)/K. 

The K-band receiver on the 70-m antenna measured only the 
left circular polarisation. The bandpass was centred mid-way be- 
tween the NH3 (1,1) and (2,2) transitions and the correlator was 
configured to have a 64 MHz wide bandpass divided into 2048 
channels. Both the NH3(1,1) and (2,2) spectra, including satel- 
lite lines, fell well inside the 850 km s usable range and were 
observed at a velocity resolution of ~ 0.4km s . 

The data were reduced in an ad-hoc pipeline built using the 
spc, gildas and miriad packages, spc was used initially to quo- 
tient and resample the spectra to a common Vlsr rest-frame. 
The spectra were inverted and frequency-shifted in the class 
package. As the spectral baselines were poor, a polynomial of 
order 5-9 (in extreme cases) was fit to the line-free channels 
before spectra from individual integrations were averaged to- 
gether. The data were then assembled into a 3D data-cube us- 
ing the cube command in the greg package and resampled onto 

2 http://www.atnf.csiro.au/computing/software/ 



a finer spatial pixel grid using the filLcube command. Finally, 
the data were output as fits cubes and miriad was used to restore 
missing header information, such as the beam-size and data-unit. 

3. Results 

Molecular emission was detected from all of the clumps identi- 
fied in the SIMBA map. The following subsections describe the 
results from the ATCA, Tidbinbilla and Mopra telescopes sepa- 
rately. 

3.1. ATCA results 

We detected NH 3 (1,1), and NH 3 (2,2) across the whole fila- 
ment. NH3 (4,4) was not detected down to a sensitivity limit of 
~ 8 mJy/beam (1.0K) when the data were smoothed to a veloc- 
ity resolution of 1 km s . Th e J,K = (4,4) level is excited by gas 
above 200 K (see Figure 1 in iHo & Townesl 19831) . hence no sub- 
stantial reservoirs of hot gas are present in NGC 3576. Figures|2} 
a and -b present the NH3 (1,1) and (2,2) integrated intensity maps 
of the region. The NH3 (1,1) map utilises data from both the H75 
and EW367 array configurations and has a resolution of ~ 1 1". 
Only H75 data exist for the NH3 (2,2) transition and the resulting 
map has a resolution of ~ 23". For reference, the 5 per cent con- 
tour (0.3 Jy/beam) of the SIMBA 1.2-mm continuum emission 
is plotted in white. The thick black line is the 3 per cent con- 
tour (0. 1 1 Jy/beam) of the 23-GHz continuum emission from the 
central Hn region. We see the NH3 emission has approximately 
the same morphology in both transitions and closely follows the 
structure of the 1 .2-mm dust emission, except immediately to the 
east of the central H11 region (SIMBA clumps S2 & S3). The 1 .2- 
mm emission in these clu mps is over 50 pe rcent contaminated 
by free-free emission (see lAndre et al.tf 2008). Sample NH3 (1,1) 
and (2,2) spectra from peak positions on the integrated intensity 
maps (labelled Sla-S5b) are presented in Figures |2]-c and -d. 
Both of the NH3 transitions exhibit the classic 'five-finger' pro- 
files at all positions. Some asymmetries between the NH3 (1,1) 
satellite lines are evident towards all positions, however, for the 
most part the spectra appear optically thin, exhibiting optical 
depths between r = 0.1 and t = 2.0. Linewidths vary between 1 
and 2.5 km s _1 , and peak brightness temperatures range between 
4.4 and 8 K for the NH 3 (1,1) line and between 1.2 and 2.6 K for 
the NH3 (2,2) lines. Individual line profiles in the western arm of 
the filament show clear evidence of blending and multiple com- 
ponents, especially towards positions S3b, S3c, S4b and S5a. 

23 GHz continuum was observed simultaneously with each 
of the NH 3 lines. To make the final map we combined data 
from the EW367, H75 and 750D arrays, including the con- 
tinuum visibilities from the line-free channels in the narrow- 
band spectral line windows. No 23 GHz continuum emission 
was detected outside of the central Hn region, down to a sen- 
sitivity limit of ~ 0.5 mJy/beam. Figure [3]presents the image of 
the Hn region at 23 GHz. 'x' and '+' symbols mark the po- 
sitions of known 6.67 G Hz methanol masers and 22 GH z wa- 
ter masers, respectively dNorris et al.lll993llCaswellll2004i) . The 
co ntours are sim i lar to the central part of the 3.4-cm map made 
by Ide Pree et alj ( 1 1999b (see Figure \Qb), with a single strong 
peak at ll A ll m 51.08*, -61 d 18 m 50.0 s (J2000) and diffuse emis- 
sion extending eastward into SIMBA clump S2. The 23 GHz flux 
density at the peak is 3.6 Jy/beam and we measure a total inte- 
grated intensity of approximately 42.6 Jy. By comparison, the 
total flux at 3.4-cm is measured t o be 71 Jy. Our obse rvations 
sample the wv-plane less well than Ide Pree et alj ( 1 1999b and we 
likely filter out a significant fraction of the extended emission. 
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Fig. 2. (a) Integrated intensity map of the NH3 (1,1) emission (thin black and white contours, and greyscale) made using the com- 
bined ATCA H75 and EW367 data. Contours are at 1.6, 7.9, 13.8, 20.0, 25.7, 31.5, 37.4, 43.6 and 50.1 K. The extent of the 23 GHz 
continuum emission from the central H11 region is ma rked by the thick black contour. The 1.2-mm dust emission is outlined by a 
single white contour (~ 3 Jy/beam, from Hil l et al.l2005l) and matches the NH3 emission closely, except in the central clump. Crosses 
(+) mark the positions of 22 GHz water masers while ' x' symbols mark the 6.67 GHz methanol maser sites. IRS sources are marked 
by circles, (b) Integrated intensity map of the NH3 (2,2) emission made using data from the ATCA H75 array, (c) Sample NH3 (1,1) 
spectra measured from the positions labelled Sla-S5b. Spectra have been offset from zero baseline for clarity. The brightness 
temperature ratios between the main and satellite lines indicate that the clumps have low to moderate optical depths, (d) Sample 
NH3 (2,2) spectra measured from the same positions. 



Our values are consistent with free-free emission from optically- 
thin ionised gas (S v °c v~ 01 ) if we are missing ~34 percent of 
the flux detected at 3.4-cm. 

We also searched for 22 GHz water maser emission across 
the filament. New masers were detected at six positions, as well 
as at the three previously known sites adjacent to IRS 1, IRS 3 
and IRS 4. Figure |4]presents a cumulative spectrum from a sin- 
gle baseline showing the majority of the maser lines detected. 
Two groupings of lines are observed, the first centred at ap- 
proximately -115kms _1 (with two components at -130 and 
-lOOkm/s) and the second centred at -30kms _1 . The intense 
lines at -130 and -100 km s -1 have peak amplitudes of ~70 



and ~ 900 Jy/beam, respectively, and derive from maser sites to- 
wards the central Hn region. Maser sites scattered across the 1 .2- 
mm continuum filament contribute to the group of lines between 
-60 and -lOkms 1 . Figure |5]presents the maser spectra and il- 
lustrates the positions of individual maser sites superimposed on 
the NH3 (1,1) zeroth-moment image. Precise coordinates, inten- 
sities and velocities are noted in Table [3] 

3.2. Tidbinbilla results 

NH3(1,1) and (2,2) integrated intensity maps made using the 
Tidbinbilla 70-m telescope are presented in Figure [6] Individual 
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Maser RA Dec Velocity Range Peak Vlsr Peak Reference. 



(J2000) (J20Q0) (kms" 1 ) (kms' 1 ) (Jy/beam) 



Ml 


11" 12"' 16.72' 


-61"17"'47.0 5 


-27 to -19 


-24.3 


5.1 


New 


M2 


11* 12'" 13.95 1 


-61 rf 17'"30.0 I 


-22 to -18 


-20.7 


3.3 


New 


M3 


ll*12 m 00.76 s 


-61 rf 18 m 31.0 s 


-41 to +15 


-25.5 


6.5 


New 


M4 


ll*ll m 53.40 s 


-61 rf 18 m 24.0 I 


-51 to +17 


-31.4 


46.5 


Caswell 2004 


M5 


ll"ll'"50.07 J 


-efis^SLO* 


-110to+12 


-101.5 


69.3 


Caswell 2004 


M6 


ll*ll m 39.80 s 


-61 rf 19 m 41.0 J 


-27 to -20 


-20.9 


7.8 


New 


M7 


11"11'"39.52 I 


-61''20'"17.0 I 


-25 to -17 


-18.4 


1.3 


New 


M8 


ll h U m 34.52 s 


-6\ d 2\ m \6.W 


-31 to -3 


-27.7 


24.8 


New 


M9 


11*11'"56.64 I 


-ei'^oi.o* 


< -134 to -121 


-129.9 


670.5 


Caswell 2004 



l 1 1 r 




-140 -120 -100 -80 -60 -40 -20 20 

Fig. 4. Cumulative H2O maser spectrum from a single baseline showing lines from across the filament. The two broad maser lines at 
-130kms _1 and -102 km s -1 derive from positions M9 and M5, respectively, towards the central Hn region (see Table[3]l. The bulk 
of the maser lines concentrated near the systemic velocity at -24kms are scattered across the filament. The units of amplitude 
are arbitrary. 



maps are centred on the same position as the ATCA fields, how- 
ever, technical difficulties meant that we were unable to finish 
one column of field S3, covering the Hn region. It is immedi- 
ately evident that the morphology of the emission is similar to 
that seen in the the ATCA NH 3 (1,1) data. 

We initially intended to use data from the 70-m antenna 
to 'fill in' the missing short-spacings in our ATCA data. For 
the merging to be successful, the emission in the single-dish 
dataset must be imaged out to its extremes. Unfortunately, there 
is still considerable flux at the edge of all fields and, due to 
time constraints and the experimental nature of the observations, 
we were unable to enlarge the maps. Instead, the Tidbinbilla 
data is used to estimate the fraction of missing flux due to ex- 
tended emission in each of the ATCA fields. To do this we 
smoothed the final ATCA maps to the same resolution and 
pixel scale as the Tidbinbilla maps and compared the integrated 
intensities measured under the same aperture. The percentage 
ATCA/Tidbinbilla integrated intensity ratios are presented in 
Table [4] Calibration error on the Tidbinbilla data is approxi- 
mately 30 percent and the measurements are consistent with less 
than 10 percent missing flux in the ATCA fields. The estimates 
assume a beam-efficie ncy of 48 per cent for the 70-m telescope 
dGreenhill et alj|200l and a (Jy/Beam)/K conversion factor of 
0.951. 

3.3. Mopra results 

Figure Q presents integrated intensity maps and sample spec- 
tra for all six molecules detected by Mopra towards NGC 3576. 
The data have been smoothed to a velocity resolution of 
0.4kms _1 and have a spatial resolution of 40". Sample spec- 



Table 4. ATCA NH 3 percentage flux detected. 



Field 


Percentage Flux 




Detected 


SI 


94 


S2 


78 


S3 


80 


S4 


91 


S5 


70 



tra presented alongside each map were extracted from the data- 
cubes at the peak positions of the SIMBA clumps (labelled SI - 
S5) identified in the 1.2-mm continuum map and marked by 
filled circles in Figure [7] 

As with NH3, the distribution of the other molecular lines 
closely follow the 1.2-mm emission except to the east of the 
central H11 region. 13 CO (1 -0) and HCO + (1-0) trace more ex- 
tended gas than the other lines, and their spectra appear moder- 
ately optically thick compared to their isotopologues, C ls O (1 - 
0) and H 13 CO + (1-0). N 2 H + and CS both trace dense gas, how- 
ever N2FF exhibits optically thin line profiles (t < 0.1) at all 
positions while CS appears optically thick in places. The rel- 
ative intensity of the emission between species varies signifi- 
cantly over the extent of the cloud, likely due to differences in 
the chemistry of the clumps. We discuss this possibility further 
in Section |531 

At position S2 the spectra of all species display three line 
components, offset in velocity by ~4kms _1 . Spectra sampled 
at the other positions appear as a single line profile. These 
lines are not well fit by single Gaussians except at position 
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Fig.5. Positions and spectra for the 22 GHz H 2 masers detected in our ATCA EW367 array data. Six masers (Ml, M2, M3, M6, 
M7 and M8) are new detections. Greyscale and black contours are the NH3 (1,1) integrated intensity map. The white line is the 
~ 3 Jy/beam contour from of the 1 .2-mm continuum map. Individual spectra are sampled from the peak maser pixel in the cleaned 
images of five ATCA fields. 
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Fig. 6. NH3 (1,1) and (2,2) integrated intensity maps made using the 70-m Tidbinbilla telescope (greyscale and contours). Contours 
levels are at 10 percent intervals. The single thick black contour shows the extent of the 1.2mm continuum emission dHill et a!l 
120051) . The grid inset into the NH3 (1,1) map illustrates the spectra which contributed to the map of the S5 field. 



S5 and are possibly composed of several blended features. 
Taking the optically thin N2H + as an example, the peak Vi sr 
of the emission in the western arm (SIMBA clumps S3-S5) 
has a mean velocity of -24.1+0.5kms _1 , while the east- 
ern arm (SIMBA clumps SI and S2) has a mean velocity of 
-25.0+0.9 km s _1 . No large velocity gradients are apparent, but 



the two arms are offset in velocity by ~ 1 km s _1 . The full-width 
half-maximum N2H + linewidths ranges between approximately 
1 and 3kms _1 across the filament, with the highest values oc- 
curring adjacent to the Hn region. 
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Fig. 7. Integrated intensity (zeroth-moment) maps of NGC 3576 in the six molecular transitions observed with the Mopra telescope 
(left panels). Sample spectra, extracted from the data cube at the positions SI -S5, are presented to the right of each map. 13 CO is 
generally optically thick and is a good tracer of the extended molecular envelope. HCO + and H 13 CO + trace similar gas densities to 
CO, while CS and N2H + trace denser gas. 



4. Derived Properties 



4.1.1. Kinetic Temperature 



In this section we derive physical properties, such as temperature 
and density, from the molecular line data. The detailed method- 
ology used in the calculations is presented in the appendix. We 
discuss the results in Section 5 where we analyse the morphol- 
ogy, kinematics, relative chemical abundances and star-forming 
content in detail. 



4.1. Physical properties from NH } 

The rotation-inversion transitions of NH3 have been used exten- 
sively in the literature to derive the kine tic temperature and col- 
umn densities of molecular clouds (e.g ., Ungere chts et al.lfl986t 
ICesaroni et aTlll992t iBoiirke et all fl 



1995I) . Due to its quantum 



mechanical properties, the NH3 inversion spectrum is split into 
multiple hyperfine components from whose brightness ratios the 
optical depth may be derived directly, hence removing a major 
assumption from the calculation of rotational temperatures. 



The modified rotational diagram method used to calculate rota- 
tional j£m£erature (T jb1 )Jto^^ desc ribed in detail 
bylHo etal.1 dl977l) . lUngerechts etalJ dl980l) and iTownes et all 
d 19831) . We provide a brief summary of the procedure in the 
appendix, specifically for the NH3(1,1) and (2,2) transitions. 
The molecular constants used in the calculations are collected 
in Table O 

Figure [8}a presents a kinetic temperature map of NGC 3576 
made using the NH 3 (1,1) and (2,2) data from the ATCA H75 
array observations only. At a resolution of ~ 23", the map re- 
flects the beam aver aged tempera ture of gas above a density 
of n=s 2 x 10 4 cirr 3 dSwaddll989t) . Higher resolution observa- 
tions would doubtlessly reveal regions with hotter or cooler tem- 
peratures than average. An obvious temperature gradient exists 
across parts of the filament. Figure (8}b plots kinetic tempera- 
ture as a function of angular offset from the peak of the ionised 
emission. The data were sampled at 12 positions, correspond- 
ing roughly to the peaks of the NH3 clumps in Figure [2] The 
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Fig. 8. (a) Map of the kinetic temperature derived from the ratio of the NH3 (1,1) and (2,2) lines. The beam size is 10.6" x 11.9". 
Kinetic temperatures range from >40K, in the centre of the map, to 10K at the extremes, (b) Kinetic temperature and (c) NH3 (1,1) 
linewidth as a function of angular offset from the peak of the Hn region. Crosses represent positions eastwards of the Hn region and 
squares positions to the west. A distinct temperature gradient is evident in the western arm of the filament (top-left). 




Right Ascension (J20O0) 

Fig. 3. 23-GHz continuum emission from the Hn region at 
the centre of the NGC 3576 complex (contours) overlaid on 
the 3.6/im SPITZER IRAC image (greyscale). The contin- 
uum image was made by combining multi-frequency data 
from the ATCA H75, 750D and EW367 arrays, including the 
line-free channels from the NH3 observations. Contours are 
at values of 0.01, 0.11 0.45, 0.85, 1.35, 1.65, 2.00, 2.48, 
2.82, and 3.25 Jy/beam. '+' and 'x' symbols mark H2O and 
CH 3OH masers, respectively , while circles mark the positions of 
the iFrogel & Perssonl (Q974) IRS sources. The total integrated 
flux density is ~42.6Jy. 



median precision on the kinetic temperature values is approxi- 
mately 2.4 K and stems from the uncertainties in the Gaussian 
fits. We find that the kinetic temperature is higher towards re- 



gions adjacent to the Hn region. In the outlying positions the 
kinetic temperature is ~ 12K, increasing to >30K towards the 
centre. Some hot spots with temperatures above 40 K exist to 
the east of the Hn region (between RAs of 11 ''12" ! 00 s and 
1 l d 12 m 15 s ), but the signal-to-noise ratio in this part of the map is 
poor, hence the derived temperatures are more uncertain (~ 8 K). 
Derived temperatures above ~ 30 K have a large associated un- 
certainty because of assumptions made du ring the calculations 
(see Appendix |A~TT>. However, below 20 K iTafalla et al.1 (2004) 
find that kinetic temperatures calculated under the same assump- 
tions are accurate to better than 5 percent. 

4.1 .2. NH 3 column density and core mass 

Assuming local thermal equilibrium (LTE), the total column 
density of NH3 may be found directly from the integrated inten- 
sity of the J,K = (1,1) line via Equation lA.5l in the appendix. An 
estimate of the total mass of gas may then be determined from 
the column density assuming a constant relative abundance to 
H 2 . 

In the ~ 11" resolution NH3 (1,1) map we begin to resolve 
substructure within the 1.2-mm clumps imaged by SIMBA. 
We have attempted to decompose the emission using the 
fellwalkeiO routine, fellwalker attempts to divide regions of 
emission into cores by searching for positive gradients in the dat- 
acube. The algorithm considers each pixel in turn above a lower 
brightness cutoff and 'walks uphill', following the steepest as- 
cent until an isolated local peak is reached. All pixels visited 
are assigned to the same core, which may already exist. Finally, 
cellular automata fill any holes in the cores and clean up the 
edges by replacing each cores' index with the the most common 
value occurring within a 3x3x3 pixel cube. The signal/noise ra- 
tio was not constant across the mosaiced NH3 (1,1) map, so the 
fellwalker routine was run on the five SIMBA regions individu- 



3 The FELLWALKER routine is part of the STARLINK package main- 
tained at http://starlink.jach.hawaii.edu/ 
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Fig. 9. (a) The spatial distribution of cores detected in the NH3 (1,1) data cube using the FELLWALKER routine overplotted on the 
integrated intensity image. Circles mark the positions of the cores and the radii are scaled to the core mass, as reported in Table [5] 
(b) Histogram showing the distribution of core masses assuming a single value of [NH3] / [H 2 ]=lxl0~ 8 . (c) Histogram showing the 
distribution of core masses corrected for the variable abundances calculated in Section l5~2l 



ally (using the mosaicked data-cube), starting at the same lower 
cutoff of 2-0", but slightly different tuning parameters. The re- 
sults were inspected by eye and compared to the original data- 
cube for consistency. Cores with fewer pixels than the area of 
the synthesised beam were omitted as unreliable detections. 

Twenty-five believable cores were found, whose positions 
are plotted in Figure|9}a. The mass of NH3 in each core was cal- 
culated from the sum of the integrated intensities of the individ- 
ual pixels in the core. We assumed excitation temperatures cor- 
responding to those derived from the 23" resolution H75 data. 
The total core mass was then e stimated assuming an abundance 
ratio [NH 3 ] / [H 2 ] of 3.0xl(T 8 dWang et alJ fcOOS). and a correc- 
tion factor of 1.38 for the abunda nce of heliu m and heavier ele- 
ments in the interstellar medium (Allenl [19731) . Table [5]presents 
the properties of the detected cores and the distribution of core 
masses is illustrated in Figure|9]-b. In Section l5T2l below, we have 
calculated the NH3 abundance at different positions in the cloud 
via a comparison to 450 //m data. Column 1 1 presents the core 
masses corrected for these results. 

Virial masses were calculated from the velocity width and 
average radius of the NH3 cores reported by fellwalker. 
Neglecting support from magnetic fields or internal heating 
sourc es, the virial mass o f a simple spherical system is given 
bv dMacLaren et al.lll988l) : 

M vir = /fcrAV 2 , (1) 

where M v ; r is the core mass in M , r is the deconvolved ra- 
dius of the cloud in parsecs, AV is the FWHM velocity width 
in kms 1 and k is an empirical constant depending on the as- 
sumed density distribution p(r). For p(r) oc l/r~ 2 , k- 126 and 
for p(r) = constant, A: = 210. Enhanced linewidths (e.g., due to 
optical depth effects or blending) will cause us to overestimate 
the virial mass and the values quoted here should be considered 
upper limits. We note that t varies between 0.1 and 2.0, so line 
broadening due to optical depth effects is negligible. Columns 12 
and 13 of Table [3] present the virial masses derived assuming 
constant and 1/r 2 density profiles, respectively. 



4.2. Physical properties from other molecules 

The molecules observed with Mopra towards NGC 3576 were 
chosen specifically to probe the physical conditions in different 
parts of the cloud. 13 CO (n crit « 6x 10 2 cirT 3 , see Table IA.lt 
and C ls O emission traces the gas in the envelope surrounding the 
dense filament. HCO + and H I3 CO + have higher critical densities 
(n C rit ~ 3 x 10 5 cm 4 ) and so probe conditions at greater depths, 
but are also proven t racers of outflows and bulk gas motions (e.g. 
iRawlings et al.l2004l) . CS and N 2 H + (n clit * 2x 10 5 ctrT 3 ) have 
comparable critical d ensities but in practise have been shown to 
trace dense gas (e.g., iPirogov et al.ll2003l lEvans|[l999l and ref- 
erences therein). N 2 H + especially, is considered to be useful in- 
dicator of the coldest (< 10 K) and most dense regi ons, where it 
is pre dicted to be the most abundant gas-phase ion dCaselh et alj 
1995). This is due to the depletion of its main destruction part- 
ner, CO, onto the dust grains. 

Assuming LTE conditions, we independently solved for the 
optical depth, excitation temperature and total column density of 
CO and HCO + via the procedure detailed in the Appendix |A. 3 1 
Figure [10] presents the final column density maps towards 
NGC 3576. The 13 CO peak optical depth ranges from 0. 1 to 4.4 
and the excitation temperature from 6.4 to 36 K. Peak HCO + op- 
tical depth ranges from 1.9 to 13.1 and excitation temperatures 
from 4.2 to 13 K. We were unable to independently estimate the 
excitation temperatures of the CS and N 2 H + transitions and so 
adopted temperatures derived from the NH3 data. This is a rea- 
sonable assumption for N 2 H + as we expect both nitrogen bearing 
compounds to be intermixed in the gas phase. We were forced to 
assume that CS was optically thin, which may not be valid as 
some of the line profiles appear saturated. 

We calculated the total mass of the molecular cloud by as- 
suming th e following abundance ratios to H 2 : [ 13 CO]/[H 2 ] = 
6.0xl0" 6 dGoldsmith et ail 1 19971 assuming an abundance ra- 
tio of [ 12 CO]/[H 2 ] =45 at a g alactocentric radius of 8 kpc 
- see lLanger & Penziasl [1990) . rHCO + l/rH 2 .l = 2.0x 10" 9 
dZhuetal.1 120071) . fCSl/TH,! = 1 . 1 x 10~ 9 dPirogov et al.l [2007h. 
and [N 2 H + ]/[H 2 ] =5.0x lO" 10 dPirogov et al.ll2007l) . The values 
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Table 5. Properties of the NH3 (1,1) cores found using the fellwalker routine in the ATCA data. 



ID 


RA 

(J2000) 


Dec 
(J2000) 


ARA" 

(") 


ADec" 
(") 


AV 
(km/s) 


Tmb 
K 


2Tmb dv 
(Kkm/s) 


N NH, 
(x 10 15 cnT 2 ) 


M LTE / 
(Mo) 


M LTE J 
(Mo) 


M viri d 
(Mq) 


— TT f — 

M vir2 T 
(Mo) 


Tkin 

K 


r 


1 


11:12:21.30 


-61:17:11.5 


12.4 


4.9 


1.3 


20.4 


2803.6 


20.8 


82 


79 


18 


11 


15 


1.5 


2 


11:12:16.70 


-61:17:39.0 


10.0 


8.9 


1.9 


23.0 


4989.6 


36.7 


144 


166 


41 


25 


20 


0.7 


3 


11:12:14.18 


-61:17:21.2 


6.4 


8.4 


2.4 


18.5 


2427.1 


18.5 


73 


168 


53 


32 


22 


0.7 


4 


11:12:08.91 


-61:18:02.9 


6.8 


6.1 


1.2 


6.1 


485.6 


3.9 


15 


15 


11 


7 


33 


0.1 


5 


11:12:05.25 


-61:18:49.3 


18.0 


8.8 


1.7 


5.6 


1538.2 


11.3 


44 


44 


48 


29 


19 


0.1 


6 


11:12:04.76 


-61:18:24.8 


7.1 


1.6 


0.8 


8.5 


131.7 


1.2 


5 


5 


3 


2 


47 


0.1 


7 


11:12:01.53 


-61:18:58.7 


5.8 


4.0 


1.0 


6.0 


232.1 


1.8 


7 


7 


6 


3 


25 


0.1 


8 


11:11:52.45 


-61:18:30.2 


4.2 


3.7 


0.6 


23.9 


185.2 


1.6 


6 


47 


2 


1 


31 


0.1 


9 


11:11:51.91 


-61:18:22.2 


7.3 


4.8 


1.0 


10.0 


668.9 


5.8 


23 


171 


8 


5 


31 


0.1 


10 


11:11:48.39 


-61:18:54.1 


2.8 


3.7 


1.2 


28.4 


809.4 


7.0 


27 


91 


6 


3 


40 


0.1 


11 


11:11:46.05 


-61:18:43.4 


8.9 


9.4 


1.7 


17.8 


3990.0 


32.5 


128 


426 


32 


19 


33 


0.1 


12 


11:11:45.66 


-61:18:52.0 


9.5 


9.3 


0.9 


9.4 


1051.4 


8.6 


34 


112 


10 


6 


35 


0.1 


13 


11:11:41.86 


-61:19:32.4 


12.8 


7.7 


1.9 


23.2 


4935.5 


37.8 


148 


62 


44 


27 


26 


0.4 


14 


11:11:39.38 


-61:19:44.2 


8.0 


5.6 


1.5 


30.2 


4757.6 


35.5 


139 


59 


20 


12 


23 


0.7 


15 


11:11:38.90 


-61:20:15.6 


3.3 


4.1 


0.8 


23.3 


647.6 


4.8 


19 


16 


3 


2 


21 


1.2 


16 


11:11:38.11 


-61:21:55.7 


3.9 


4.6 


1.6 


4.3 


222.1 


1.8 


7 


14 


13 


8 


- 


0.1 


17 


11:11:37.94 


-61:21:06.5 


1.3 


2.8 


0.8 


85.9 


595.2 


4.5 


18 


35 


1 


1 


15 


1.2 


18 


11:11:37.92 


-61:19:35.1 


6.6 


4.4 


1.0 


17.5 


975.2 


7.3 


29 


12 


7 


4 


24 


1.0 


19 


11:11:36.49 


-61:20:07.0 


7.6 


12.3 


1.7 


25.9 


6502.5 


47.9 


188 


166 


36 


22 


17 


1.0 


20 


11:11:35.62 


-61:21:30.0 


7.9 


2.6 


0.6 


20.5 


422.7 


3.2 


12 


25 


3 


2 


15 


2.2 


21 


11:11:33.11 


-61:21:11.6 


13.2 


8.4 


2.4 


27.9 


8117.8 


65.7 


258 


516 


78 


47 


15 


1.9 


22 


11:11:32.84 


-61:19:53.0 


1.2 


0.3 


0.8 


321.3 


231.1 


1.7 


7 


6 


1 





17 


1.0 


23 


11:11:27.13 


-61:21:07.3 


5.6 


2.0 


0.7 


32.3 


517.0 


3.9 


15 


30 


2 


1 


16 


1.6 


24 


11:11:24.25 


-61:21:18.7 


5.2 


6.2 


0.8 


4.4 


109.9 


0.9 


3 


7 


5 


3 




0.1 


25 


11:11:24.03 


-61:21:25.3 


6.5 


6.6 


0.9 


8.6 


534.0 


4.2 


17 


33 


7 


4 




0.1 



" Deconvolved full- width half-maximum (FWHM) core size asumming an average beam FWHM of 1 1.2". 

13 LTE-mass derived from the integrated intensity of NH 3 (1,1) assuming a relative abundance [NH 3 ]/ [H 2 ] = 3 x 10~ 8 . 

y Corrected LTE-mass, derived using the [NH3]/[H2] abundance ratios calculated in Section l5T2l via a comparison between NH3 and 

450 yt/m p- ArTeMiS data. 

6 M v ii derived assuming a constant density profile. 
1 M vir derived assuming a 1/r 2 density profile. 



for total cloud mass derived from the CO, HCO + and N2H + ob- 
servations are all approximately 8,300 M . The mass derived 
from CS is 1 2,700 M . These values assume calibration onto 
the Mopra 'extended beam' temperature scale, which includes 
power received from the first sidelobe and is appropriate for data 
more extended than 80" (see Section 2.1 of lLadd et al. 20051). 

5. Analysis & Discussion 

5.1. H 2 masers 

The three well-known water maser sites within the central Hn re- 
gion (labelled M4, M 5 and M9 in Figure [5]l have been in- 
vestigated by ICaswelll (I2004I) . who summarises their proper- 
ties. At the M4 site (G291. 274-0.709) the H 2 maser coex- 
ists with 6.67 GHz CH 3 OH and 1.67 GHz OH masers. In the 
H2O maser spectrum we detect bright maser features rang- 
ing from -37kms _1 to 23kms~ 1 , in agreement with Caswell, 
however, we also detect weak (<5Jy) lines out to velocities 
of -Slkms- 1 and +20kms- 1 . The M5 site (G29 1.270-0.7 19) 
hosts a weak CH3OH maser at -26kms~' and a bright 
H2O maser a t a ve locity of -102kms~ 1 . First reported by 
ICaswell et all dl989l) . we find this H2O maser has an unusu- 
ally broad linewidth (2.9kms _1 ) and a peak flux density of 
69 Jy, consistent with the original observations. The Vlsr of 
this spectral feature is reported to have become increasingly 
negative since its discovery in 1989 when it had a velocity of 
-88kms _1 . This has been interpreted as being due to acceler- 
ation of the emitting gas, perhaps in a high velocity outflow 



(ICaswelll2004t ICaswell & Phillipsl2008l) . Alternatively, we sug- 
gest that the -102 kms -1 line may be a new maser feature not 
associated with the component seen at -88kms~', which may 
have been quenched in the intervening time. The final known 
site, M9 (G29 1.284-0.7 16), exhibits a single intense maser at a 
velocity of -130 km s _1 and with a linewidth of 4.8kms~ I . We 
measur e an intensity of 670 Jy in contrast to the earlier value of 
945 Jy (ICaswell et al.l 1989b . perhaps reflecting a real decrease in 
the intensity over the intervening time. We note that this maser 
feature falls on the edge of our bandpass, which may make the 
flux density measurement unreliable. 



The six new maser sites reported here (Ml, M2, M3, M6, 
M7 and M8) are located in the 'arms' of the filament. As can 
be seen from Figure [5] all sites are within a few arc-seconds of 
NH3 emission peaks. With the exception of M8, it is notable that 
the spectra of these new masers exhibit, at most, two bright fea- 
tures. These features have velocities close to the systemic veloc- 
ity and intensities below 10 Jy. In contrast, the maser spectrum of 
M8 has four strong peaks spread over ~ 35 km s _1 , the brightest 
of which has an intensity of 25 Jy. The mid-infrared image of the 
host SIMBA clump (S5) also exhibits some nebulous emission, 
on the south-west edge (see Figure QJa), while the dense molec- 
ular gas tracers (see Sections 13.11 13.21 and 13. 3t show a 'notch' 
in the emission at the same site. We interpret these features as 
evidence of a deeply embedded young stellar cluster evacuating 
an open-ended cavity. 
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Fig. 10. Column density maps of the four species targeted by the Mopra telescope. Contours are set at 10 per cent intervals, beginning 
at 10 per cent, except in the HCO + column density map, where they start at 30 per cent of the peak value. 



5.2. Core masses 

Core masses derived from NH3 rely on the assumption of a 
constant [NHaJ/fHa] abundance ratio, whose value may vary 
across the cloud. Independent estimates of the gas + dust-mass in 
NGC 3576 have been made bv lAndre et al.l d2008l) . who mapped 
the 450 /im thermal emission using the P-ArTeMiSQ bolometer 
camera on With a beam size of 10" the 450//m map 
has a comparable resolution to our NH3 (1,1) observations and 
the morphology of the 450 fim emission corresponds almost ex- 
actly to the integrated molecular emission. 

We expanded the analysis of lAndre et alj ([2008) by decom- 
posing the 2-D 450 yum emission into cores using fellwalker and 
deriving their individual masses via the relation: 

M = S45 °^ mD2 , (2) 

#450 Bv(Td)450 

where S450 is the measured flux density, K450 is the dust opac- 
ity per unit (gas + dust) mass column density at 450 fim, B v (Td) 
is the Planck function B V (T) at a dust temperature and D is 
the distance to the source. We set K450 = 0.04cm 2 g _1 , a value 
generally adopted in sub-mm studies of pre-s tellar cores (e.g. 
iMotte et al.ll 19981 IWard-Thompson et al.ll 19991) and appropriate 
in r egions of moderately high gas densities (around 10 5 ctrr 3 - 
c.f. iHenning et al.ll 1995b . For consistency, we set Tj to the ki- 

4 See http://irfu.cea.fr/Sap/en/Phocea/Vie_des_labos/Ast/ast_technique. 

5 Atacama Pathfinder Experiment http://www.apex-telescope.org/ 



netic temperatures derived from our NH3 data. A direct compar- 
ison between the masses derived from 450 fim and NH3 emission 
integrated over the same apertures results in estimates of the av- 
erage [NH3]/[H2] abundance ratio in each core, assuming no 
depletion. 

Eight equivalent cores were detected in the 450/mi P- 
ArteMiS and ATCA NH3 (1,1) integrated intensity maps, whose 
positions and masses are presented in Table [6] Mass values 
from 450jt/m data range from ~89M G to ~712M . Some 
450 jum cores encompass two or more NH3(1,1) cores and in 
such cases we summed the masses of the components quoted in 
column 10 of Table [5] Relative [Npy/fpy abundance values 
range from 0.4x1 0" 8 to 7.1xl0~ 8 and are distinctly different in 
adjacent cores. We note that the error on the abundance values is 
approximately a factor of 3, including the calibration uncertainty 
of the P-ArteMiS data, assumptions about the dust emissivity 
and opacity, dust-to-gas ratio, temperature and relative angular 
size. The abundance differences may also be attributed to en- 
vironmental factors. Cores S3-a and -b exhibit the lowest abun- 
dance values, but lie adjacent to the Hn region. Here they may be 
exposed to the ultraviolet-radiation field from the central cluster, 
leading to enhanced destruction of large molecules. Core S4- 
a, and to a much lesser extent S4-b, have elevated abundances 
of NH3. Interestingly, this SIMBA clump also has an enhanced 
abundance of N2H + , possibly due to depletion of CO onto the 
dust grain surface (see Section 15.31 below). In light of the vary- 
plipgidtlaHiftd_aftSi values we have attempted to correct the masses 
derived for the cores found in the 3-D NH3 data-cube. Where an 
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Table 6. Comparison of masses derived from p-ArTeMiS 450//m and ATCA NH3 (1,1) emission. 



Clump 




RA 


Dec 


Mass 


Mass 


[NH 3 ]/[H 2 f 


Average 


Aliases 


ID 


(J2000) 


(J2000) 


450;um 


NH 3 LTE 


Abundance 


Tdust 


(Andre et al. 2008) 










(Mo) 


(M Q ) 








Sl-a 


11 


12:21.09 


-61:17:09.50 


89 


91 


3.1 x 10~ 8 


15 




Sl-b 


11 


12:17.21 


-61:17:43.50 


161 


140 


2.6 x 10~ 8 


21 


Sl-Ml 


Sl-c 


11 


12:14.43 


-61:17:17.50 


193 


82 


1.3 x 10- 8 


23 


S1-M2 


S3-a 


11 


11:53.05 


-61:18:31.37 


357 


51 


0.4 x 10- 8 


40 


S3-M4 


S3-b 


11 


11:45.54 


-61:18:45.27 


682 


211 


0.9 x 10- 8 


38 


S3-C3 


S4-a 


11 


11:39.69 


-61:19:43.18 


161 


383 


7.1 x 10- 8 


24 


S4-M6 


S4-b 


11 


11:36.91 


-61:20:03.13 


247 


279 


3.4 x 10~ 8 


19 




S5 


11 


11:34.66 


-61:21:11.09 


712 


339 


1.4 x 10 _s 


15 


S5-M8 



Relative abundance of NH 3 to H 2 via a direct comparison of clump masses derived from NH 3 and 450 fim emission. 
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Fig. 11. Distributions of M correct ed/M v i r for the 25 cores found in 
the NH3 data. The M conectet j values have been derived from the 
Mlte masses scaled to the 450yum core masses. M v ; r was calcu- 
lated assuming a assuming a r~ 2 density profile. For a constant 
density profile virial masses would be 1 .7 times larger. All cores 
have ratios M correcte d / M v j r > 1 and are at least consistent with 
being gravitationally bound. Cores with ratios much greater than 
one are likely to be gravitationally unstable and may be under- 
going collapse unless supported by magnetic fields on the order 
of milli-Gauss 



equivalent 450 /mi c ore does not exist we revert to the canonical 
value of 3.0X10 -8 dWang et al .1120081) . These 'corrected' values 
are presented in column 1 1 of Table [5] 

FigureQT]compares the corrected LTE and virial masses for 
the NH3 cores in Table [5] The factor f = M correct ed/M v i ri mea- 
sures the ratio of gravitational to kinetic energy. Values less than 
one indicate that the core may be a transient structure or con- 
fined by an over-pressured external medium, while a value of 
f ~1 suggests that the core is close to gravitational equilibrium. 
Cores with f »1 are likely to be gravitationally unstable, mean- 
ing that internal magnetic fields are required to counteract the 
effects of self-gravity. We find that in all cases f > 1 (within er- 
rors) and seventeen cores have f-values between three and thirty. 
The magnetic fiel d strength necess ary to support each core may 
be calculated after lBot et al] J2007) via 



1 - ^ 



10 



10\ GMVo 



9f R 4 7T 



(3) 



where Bo = 0.3 nT is the am bient magnetic field strength in the 
Milky Way dHan et al. 2006), R is the core radius in metres and 
M is the core mass in kilograms. G and yUo take their usual val- 
ues in SI units. We calculate that field strengths ranging from 



~1 mG to ~40mG, with a median of ~4mG, are required to bal- 
ance sel f-gravity in addition to the turb ulent pressure. By com- 
parison, |Curran&^Chr^sostomou| (|2007|) conducted polarimetry 
measurements of the magnetic field strengths in fourteen mas- 
sive star forming regions and found values from < 0. 1 mG to 
~6 mG. Our values are higher on average, suggesting that some 
of our NH3 cores may indeed be collapsing. In Section |5~4l below 
we present evidence of collapse in the SIMBA clump S4. 

lAndre et"ail (|2008) estimated the envelope mass of the dom- 
inant protostellar object in each 450//m core by scaling the mass 
corresponding to the peak flux density to a uniform diameter 
of 6,000 AU, assuming a density profile p(r) ocr~ 2 . The values 
range from 21 to 45 M . They also estimated the bolometric lu- 
minosities by fitting the SIMBA, P-ArTeMiS and MSX0 data 
with a grid of s p ectral energy distribution models computed by 
iRobitaille et all d2006l 120071). When plotted on a M env versus 
Lboi diagram (see Figure 4 of lAndre et al.l l2008) the objects are 
bracketed by evolutionary tracks corresponding to final stellar 
masses of 15 to 50 M Q , implying that the cores are undergo- 
ing high-mass star-formation. The relative positions of the cores 
on the diagram also indicates their relative evolutionary ages. 
Interestingly, the cores furthest from the Hn region, at positions 
Sl-b, Sl-c, S4-a, and S5, are significantly less evolved than the 
adjacent cores at S3-a and S3-b. This evolutionary gradient is 
consistent with the hypothesis that the central Hn region has trig- 
gered sequential star-formation in the filament. 

5.3. Chemistry 

Once the physical conditions of temperature, density and 
optical-depth have been determined, we can begin to examine 
the chemistry in NGC 3576. Ratios of column density provide 
beam-averaged measures of the relative abundances between 
two species. For key molecules this ratio is tied directly to the 
physical conditions and hence to the embedded star formation 
activity. 

Figure Q~2] presents maps of the relative abundance of 
[N 2 H + ] / [CO] and [N 2 H + ] / [CS] made by taking the ratios of the 
column density maps. The abundance ratios are beam-averaged 
and assumes all of the gas is above the critical density for each 
species. Immediately clump S4 stands out in the [N2H + ]/[CO] 
map, as it has an over-abundance of N 2 H + compared to the 
other SIMBA clumps. The difference is not as pronounced in 
the [N 2 H + ] / [CS] map, however, we note the optical depth of 
CS is not known and hence there are significant uncertainties 



6 Midcourse Space Experiment (MSX) 
http://irsa.ipac.caltech.edu/applications/MSX/ 
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Fig. 12. Maps of the relative abundance of [N2H + ] / [CO] (left) and [N2H + ] / [CS] (right) made from individual column density 
maps. Greyscale and contours are set at 10 per cent intervals, starting at 30 per cent of peak. The N2H + abundance is significantly 
enhanced towards clump S4. This difference is especially pronounced in the [N2H + ] / [CO] map but is not as clear when comparing 
N2H + and CS. We note that the column density of CS is the least well known of the three species. 



in the column density across the map. A comparison between 
the N2H + and HCO + column densities also reveals the same en- 
hancement in clump S4, which we believe to be a real chemical 
difference. Because of the low signal-to-noise we do not show 
this map. 

The abundance of N2H + is predicted to be enhanced in cold 
(< 10 K) or dense clumps. This is because its parent molecule 
N2 is one of the least affected by the condensation process and 
N2H + survives in the gas phase, at least for densities in the range 
10 5 — 10 6 cm -3 . In such environments its main destruction part- 
ner CO depletes ont o the grains and the abundance of N2H + con- 
sequently increases (Bergin et al. 2002). Interestingly, 13 CO and 
C 18 emission are noticeably absent from clump S4, providing 
further confidence in this interpretation. Conditions such as these 
are found in starless cores, but particularly in cold collapsing 
clumps evolving towards forming stars. Clump S4 has an av- 
erage kinetic temperature of ~ 25 K according to our analysis 
of NH3 data. This seems high for an infrared dark cloud (see 
Figure [T}a) undergoing collapse, however, CO depletion may 
still occur at these temperatures if the density is high enough 
dAikawa et al.ll200"lh . In the following sub-section we will ex- 
amine the evidence for inflowing gas motions in this clump. 

5.4. Morphology and kinematics 

We have examined the molecular emission for evidence of veloc- 
ity gradients indicating bulk gas motions, expansion or contrac- 
tion. Figure [13] presents position-velocity (PV) diagrams made 
using the 13 CO Mopra data and annotated on the inset inte- 
grated intensity map. Panel 1 was made using the kshell tooQ and 
shows the azimuthly averaged intensity as a function of veloc- 
ity and angular offset from position 11 /, 11'"58.0 5 -61 d 18 m 35.0 s 
J2000. The 'C shaped profile is a classic indicator of an ex- 
panding shell of gas, in this case centred approximately on the 
Hn region. Towards the centre position (zero offset) the spa- 
tial pixels sample gas moving directly towards and away from 
the observer, hence the average spectrum peaks at the systemic 
velocity + the shell expansion velocity. At offsets approaching 
the shell radius the pixels sample gas moving perpendicular to 



7 kshell is part of the KARMA suite of analysis software available 
from the ATNF at http://www.atnf.csiro.au/computing/software/ 



the observer, hence the emission peaks at the systemic velocity 
(-24km s _1 ). The clumpy shell seen here in 13 CO is approxi- 
mately 1 arcminute in diameter and encompasses the free-free 
emission from the H11 region. We see exactly the same C-shape 
in C ls O, although the plot is noisier. 

Panels 2-4 present the PV-slices indicated on the in- 
set 13 CO moment map (top-left). The first cut, (centred at 
11*11 B, 09.90* ) -61 rf 15 m 45.6* J2000 and at a position angle 
= 244.5°) is a longitudinal slice through the eastern arm of 
the filament and reveals a forked cross-section: two tenuous ten- 
tacles of emission merge into the eastern-most SIMBA clump 
(SI). PV-slices through the dense gas tracers (NH3 and CS) show 
the same structure, meaning that this morphology is not due to 
optical depth or chemical effects. Ionised gas from the Hn region 
overlays this part of the eastern arm, likely inter-penetrating the 
molecular emission. Combined with the existence of an expand- 
ing shell and a temperature gradient, we conclude that the Hn re- 
gion is expanding eastwards, and has played a significant role in 
sculpting these features. 

The second cut (1 l''l l'"48.60 ! , -61 i/ 20 m 57.6 s J2000, 
= 174°) follows the line of emission evident in Figure [T3l at 
— 20kms _1 . It is clear from the PV-diagram that two velocity 
components exist, separated by ~4kms~ 1 . Both appear to merge 
into the strong emission immediately west of the Hn region. 
These two velocity features are clearly visible in the third cut 
(1 1 A 1 l ffl 39.50 s , -6\ h 2\ m A6.9 s J2000, = 268°), which presents 
the velocity structure of the western arm in cross-section. The 
bulk of the emission in this arm derives from the component at 
-24km s" 1 , while the component at -20km s is not seen in 
the dense gas tracers NH3, N2H + or CS. Considering the dense 
gas only, we see the peak Vi sr of the filament is approximately 
constant between -24 to -25 km s _1 over its length. No large op- 
posing velocity gradients or discontinuities are evident and small 
variations may be explained by optical depth effects or motions 
on small scales (< 30")- We see no evidence for large scale flows 
of gas, at least in the 40" resolution Mopra data, and conclude 
that the majority of the emission derives from the same cloud 
and is not merely a projection or line-of-sight effect. The veloc- 
ity structure in the 11" resolution NH3 data is complex within 
the bounds of each SIMBA clump and will be analysed in detail 
in a separate paper. 
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Fig. 13. Position-velocity cuts through the 13 CO data-cube. The top-left panel is an integrated intensity map of 13 CO showing the 
orientations and positions of the PV-slices presented in panels 1-4. Centre positions of the slices are marked with black and white 
filled circles. Panel 1 shows the azimuthly averaged intensity as a function of velocity and angular offset from position 1 l"ll m 58.0 i 
— 61 d 18'"35.0 s J2000. The 'C shaped profile is indicative of an expanding shell of emitting gas. Panels 2, 3 and 4 present PV-slices 
along the axes of emission and in cross-section through the western arm, respectively. Several distinct filamentary structures are 
evident, merging into the bright emission adjacent to the Hn region. 



We are interested in probing for bulk gas motions in clump 
S4, where we find an over-abundance of N2H + , common in col- 
lapsing clumps. Asymmetric spectra towards star forming re- 
gions are often interpreted as indicators of inward or outward 
motions. In the case of NH3 (1,1) asymmetries may be pro duced 
by a combination of non-LTE effects and bulk motions. iParkl 
(2001) has modelled NH3(1,1) spectra over a wide range of 
physical conditions, predicting that asymmetries between the in- 
ner satellite lines indicates inward or outward flows of gas. On 
the other hand, asymmetries between the outer satellites are in- 
dicative of selective radiative trapping and non-LTE excitation 
dStutzki & Winnewisser|[l985t) . The amount by which both ef- 
fects skew the profiles increases with optical depth. Figure \14[ 
a presents a map showing the difference between the blue mi- 
nus red peak intensities of the inner satellite NH3 (1,1) lines, for 
SIMBA clump S4. Only pixels with absolute differences greater 
than 0.2 K (~ 3-cr above the spectral noise) have been plotted. 
Asymmetries on the order of ~ 10 percent are observed over 
distinct regions within the SIMBA clump. We note that there 
are few pixels with significant red-skewed spectra (values less 
than -0.2 K in Figure [l4"la) and the map is dominated by re- 
gions of blue-skewed spectra which cover angular areas greater 
than the beam. Sample spectra from the three largest regions are 



presented in Figure [T4lb. The optical depth of NH3 (1,1) in this 
clump ranges between 0.4 < r < 1.0 and we do not expect 
large asymmetries under these conditions, hence, we believe the 
differences are significant and tentatively imply inward gas mo- 
tions. 



5.5. H11 region - cloud interaction 

It is clear from the temperature map presented in Figure [8] that 
the central star-forming complex is heating the gas in its imme- 
diate surroundings. The heating is gradual on the western side of 
the filament and the temperature decreases linearly towards the 
extremes. The gradient in temperature strongly suggests that the 
Hn region is embedded within the filament and is not merely a 
line of sight projection. The temperature fluctuates between ~ 15 
and ~ 30 K in the S2 clump, immediately east of the Hn region's 
peak. Diffuse free-free emission is observed to extend into this 
region and the molecular gas has a clumpy distribution and is 
less evenly heated. It is interesting to note that there is a weak 
linewidth gradient in the western arm also, as plotted in Figure[8}- 
c . This may reflect the turbulent energy injected into the cloud 
via shocks driven by the expanding ionised gas in the Hn re- 
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Fig. 14. (a) Map of the difference between the peak intensities of the F J <-> F J" (blue) minus F [ «-» F j (re d) wwer satellite line s of 
NH3 (1,1) for SIMBA clump S4, where F is the total angular momentum quantum number (see Figure 1 in lRvdbeck et al.|[l977l) . 
Only spectra with differences greater than 0.2 K (~ 3-cr) are plotted as colourscale, while contours represent the integrated intensity map. The 
inset histogram shows the distribution of difference values over the whole of the clump and is clearly dominated by significant numbers of 
blue-skewed spectra, (b) Sample spectra from the positions marked i, ii and iii on the difference map. 
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Fig. 15. The distribution of st ars with infrared excess deter- 
mined by Mae rcker et al.l d2006l) . plotted over the molecular gas 
traced by 13 CO (1 -0) emission. 



gion. Clear evidence for this continuing expansion is seen in the 
13 CO data as an expanding shell of molecular gas surrounding 
the free-free emission. 

5.6. Infrared excess stars 

Previous studies have focused on the infrared properties of em- 
bedded stars in NGC 3576. Emission in the K (2.2yum) and L 
(3.6/mi) bands, which cannot be accounted for by interstellar 
reddening, has been inte rpreted as the signature of thermal emis- 
sion from a dusty disk. Maercke r et alj J2006) imaged the re- 
gion in J (1.3//m), H (1.7//m), K and L, and found 113 stars 
with infrared excess. Their distribution is plotted in Figure [15] 
over the 13 CO (1 -0) integrated intensity map. Foreground stars 
have been filtered from the list. A dense cluster of IR-excess 



stars is centred on the Giant HII region, as reported in work by 
iPersi et ail d 19941) using Las Campanas Observatory. It is clear 
that the distribution of IR-excess stars is anti-correlated with the 
molecular gas in the dense arms of the filament. The high extinc- 
tion in this area likely makes it impossible to detect even bright 
IR-excess sources behind or embedded in the filament. This sug- 
gests that the clumps are not associated with star-formation at an 
advanced stage of development. If young stars are present in the 
filament, they are probably at such a young age that they are still 
deeply embedded and not yet visible in the near-infrared. 

5. 7. Triggered star formation in NGC 3576? 

Our observations provide an interesting glimpse into star forma- 
tion activity within NGC 3576. It is clear that the Hn region is 
embedded within the dusty filament and is exerting a powerful 
influence on its immediate surroundings. We have investigated 
whether the Hn region expansion into the ambient molecular 
cloud is consistent with triggering the formation of high-mass 
stars along the dusty filament. 

By comparing the morphology of the Hii region (traced 
by the 23 GHz continuum emission) to the dense molecular 
gas (traced by NH3, CS and N2I-F) and thermal dust emis- 
sion (traced by 1.2-mm and 450yum continuum), we see that the 
young star cluster at the heart of the complex has destroyed or 
dispersed much of the gas in the central clump. There is also 
clear evidence for an expanding shell in the 13 CO and C ls O data, 
centred on the Hn region. 

To the east of the Hn region, molecular gas is sparse and does 
not follow the morphology of the dust, while at 3.4-cm wave- 
lengths the ionised gas is observed to extend throughout most 
of field S2 (see FigureQJc). In contrast, the western edge of the 
free-free continuum emission is observed to have a well defined 
boundary. We find from our molecular-line observations that this 
edge coincides with a dense molecular clump, implying that the 
Hn region is pressure-confined on this side. At the same time, 
the Hn region appears to be heating the arms of the filament and 
the western arm exhibits a clear temperature gradient, peaking 
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Fig. 16. Schematic of the environment of NGC 3576 (right panel). 13 CO integrated intensity contours tracing the extended molecular 
cloud, overlaid on a H-a image (left panel), which highlights the large scale structure of the interstellar medium around the star- 
forming complex. The molecular emission is correlated exactly with the extincted regions of the H-a. Loops of H-a emission are 
seen to extend north and south from the vicinity of the ionised Hn region (thick black contour), likely shaped by the energetic young 
cluster at the heart of the complex. The ionised gas is expanding into the molecular cloud, sweeping up a clumpy shell of gas and 
heating the filamentary cloud in which it is embedded. Dense cores containing tens to hundreds of solar masses have formed along 
the length of the cloud and exhibit 22 GHz water maser emission, signifying star-formation in progress. 



adjacent to the Hn region. In the eastern arm, the clumpy nature 
of the gas gives rise to 'hot spots' with higher temperatures, con- 
sistent with our 'dispersal' interpretation. The gradients in tem- 
perature strongly suggest that the Hn region is embedded within 
the filament. Some of the structures visible in the eastern arm are 
reminiscent of the 'elephant's trunks' pillars in the Eagle Nebula 
(IWhiteet al.ll 19991) . 

Besides the heating evidence, there is also evidence for 
kinematic al interaction between the Hn region and the ambient 
molecular cloud. We do see evidence that the Hn region is ex- 
panding into the eastern arm, sweeping up a clumpy shell of gas. 
Such an expansion is expected in the collect and collapse sce- 
nario for triggered star formation (see Section[T|). Multiple veloc- 
ity components have been identified within the arms of the fila- 
ment (see Figure[T3Tl. Velocity gradients and dislocations may be 
the archaeological remains of clo ud-cloud collisions or quickly 
collapsing elongated clumps (c.f. Pere tto et al.ll200 6). However, 
when considering only the dense gas we see that the velocity 
structure of the gas is smooth, confined to a narrow range and 
does not support either interpretation, at least at the 40"scale of 
the Mopra data. 

Different diagnostics have been employed to identify sig- 
natures of early phases of star formation in the filaments and 
gas clumps. Water masers have alway s been associated with 
star formation (see iBeuther et alj 120021 for a summary). That 
we detect water-maser emission towards the arms of the fila- 
ment confirms the presence of shocks and hence that stars are 
forming within the molecular gas. However, unlike 6.67 GHz 
CH3OH masers, 22 GHz water masers are not a clear confir- 
mation of high-mass star formation as they are also found to- 



wards sites where intermediate mass (4-8M ) stars are form- 
ing. Water masers are thought to originate in the shocked gas 
associated with outflow s and are collisionally pumped (e.g., 
iKvlafis & Pavla kis 1992). Core masses derived from NH3 range 
from 5 M to ~ 5 16 M , with a median value of ~ 45 M . Seven 
cores have masses greater tha n 100 M and for a t ypical star for- 
mation efficiency of 0.3 (c.f. lLada & Ladall2003l) are likely to 
form stars with masses greater than 8 M . 

We have confirmed that protostellar objects are likely form- 
ing in the arms of the filament and that the Hn region is clearly 
influencing the conditions in the bulk of the cloud, as evidenced 
by the temperature gradients and gas-morphology. Whether the 
formation of the new stars is triggered by the expansion of the 
Hn region is difficult to prove categorically, but to be consistent 
requir es evidence of an 'age gradient', as argued b y iPersi et al. 
(1994). Recent 450-//m continuum observations dAndre et alj 
2008) partially resolved the 1.2-mm continuum SIMBA clumps 
and detected continuum emission counterparts to the NH3 cores. 
Based on a 'protostellar mass envelope vs. bolometric luminos- 
ity diagram' (hereafter M env -L[, \ diagram). lAndre et ail (2008) 
estimated the evolutionary stage of these candidate massive pro- 
tostars and positioned them on evolutionary tracks for high mass 
star formation. Interestingly, S3-a and S3-b, the closest candi- 
date protostars to the Hn region, appear to be the most evolved 
objects on the M em , - Lhoi diagram. In contrast, the other candi- 
date protostars, those in the arms, appear to be less evolved on 
the diagram and are possibly class 0-type protostars. The non 
detection of ionised gas toward these protostars confirms their 
very early evolutionary stage. All detected protostars will evolve 
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into stars of masses between 15 and 50 M according to the evo- 
lutionary tracks on the M env - L\,„\ diagram. 

Another signature of early phases of star formation is the 
evidence for collapsing cores. On examination of the relative 
molecular abundances, SIMBA clump S4 stands out as having an 
over-abundance of N2H + and an under-abundance of CO. Such a 
chemical signature is typical of dense and cold regions, which is 
often found towards collapsing cores. We searched for signatures 
of inflowing motions in the NH3 line profiles. Asymmetries on 
the order of ~ 10 per cent are observed between the inner satel- 
lites for a region running lengthwise along the clump S4. 

High-mass star formation is clearly underway in the dusty 
arms of the NGC 3576 filament. Evidence for interaction 
through temperature gradients as well as an evolutionary se- 
quence for the embedded protostars from the Hn region to the 
arms, is consistent with the hypothesis that the Hn region is trig- 
gering star formation in NGC 3576. The geometry of the Hn re- 
gion and associated clumps does not resemble the typical collect 
and collapse scenario, hence the observations are also in agree- 
ment with spontaneous star formation, at least in the clumps fur- 
ther away from the ionisation front. However, it is also possible 
that the Hn region is responsible for initiating star formation in 
these clumps via shock driven globule squeezing. Further inves- 
tigation is necessary to distinguish between these two cases. 

6. Summary & Conclusions 

We have mapped the molecular environment of the giant Hn re- 
gion NGC 3576 in lines of NH 3 , CO, HCO + , CS and N 2 H + . 
In addition, we have searched for UCHn regions via 23 GHz 
free-free emission and for 22 GHz water masers. FigurefTBlsum- 
marises our observational results as a schematic, overlaid on an 
amateur H-or image (J. Schedler, 2009) which illustrates the large 
scale structure of the ISM in the region. Our main conclusions 
are as follows: 

1 . We have detected molecular emi ssion from all clu mps identi- 
fied in the 1.2-mm SIMBA map dHill et al.ll2005l) . Dense gas 
tracers (NH3, N 2 H + and CS) follow the morphology of the 
thermal dust emission, except within the Hn region, where 
the 1.2-mm emission is highly contaminated by free-free 
continuum emission. Molecules with low and intermediate 
critical densities or with high abundances (CO, HCO + ) trace 
the extended envelope well. 

2. We searched for UCHn regions in the arms of the filament via 
their 23 GHz emission, but found none down to a detection 
limit of 0.5 mJy/beam. The central Hn region was mapped 
and ex hibited a similar m orphology to the earlier 3.4-cm im- 
age bv lde Pree et al.1 d!999t) . 

3. Six new sites of 22 GHz water maser emission were detected, 
all of which lie adjacent to discrete NH3 cores in the arms of 
the filament. Most of the new masers have flux densities less 
than ten Janskys, compared to the three known sites in the 
Hn region, which have flux densities between 40 - 640 Jy. Of 
all the new masers, the site in clump S5 has the greatest range 
of velocities and warrants further investigation. 

4. NH3 (1,1) and (2,2) were used to derive the kinetic tempera- 
ture of the dense gas. We find a temperature gradient exists 
in the western arm, with temperatures ranging from ~ 1 1 K 
at the extremes to > 30 K adjacent to the Hn region. Clearly, 
the Hn region is responsible for the heating. 

5. There is clear evidence in 13 CO and C 18 for an expand- 
ing swept-up shell of diameter ~ 1' completely encompass- 
ing the Hn region. 



6. The ~ 11" resolution NH3(1,1) emission was decomposed 
into twenty-five individual cores, ranging in mass from 
5M Q to 516M , when scaled to the 2-D clump masses de- 
rived from the 450/im P-ArTeMiS map. 

7. M correctei i / M v i r >1 for all NH3 cores, implying they 
are at least gravitationally bound. Seveneen cores have 
M conec ted / M v i r > 3 and to avoid collapsing under self-gravity 
require support from internal magnetic fields with strengths 
between ~1 mG and ~40mG These values are on average 
higher than measured in recent observations of massive star- 
formation regions and suggest that some of the cores may be 
collapsing. 

8. Clump S4 in Figure [TJ is over-abundant in N 2 H + and 
shows evidence for depletion of CO, consistent with the 
expected chemistry in a cold, dense clump. We examined 
our NH3 data for line-profile asymmetries indicating bulk 
gas motions and found distinct regions with blue profiles, 
tentatively indicating inflowing gas. 

9. The filament displays a complex velocity structure, with 
two main components in each arm seen in the I3 CO data. 
Molecules tracing high density gas show that a single com- 
ponent at Vlsr = -24 km s _1 dominates the cloud. 

10. The positions of stars with infrared-excess, indicative of cir- 
cumstellar disks, are anti-correlated with the dense molec- 
ular gas, implying star formation has not yet advanced far 
enough to be visible in the near-infrared. Alternatively, there 
is likely sufficient extinction to hide such stars embedded 
within the molecular cloud. 

From the above evidence it is clear that star formation is un- 
derway across the whole of the NGC 3576 filament and that the 
Hn region is influencing the physical conditions in the bulk of the 
molecular gas. Very young massive protostellar objects, some 
equivalent to low-mass class-0 protostars, are forming within 
the arms of the filament. Masses and luminosities of cores with 
1.2-mm, 450//m a nd infrared data have been characterised by 
lAndre et"ail (|2008). Objects adjacent to the Hn region appear to 
be more evolved, implying an evolutionary gradient. The Hn re- 
gion is expanding into the cloud, dissipating the gas on its east- 
ern edge and is possibly the trigger for star-formation in the near- 
est clumps. However, the observational data is also consistent 
with spontaneous star formation, at least in the outlying regions. 
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Appendix A: 
lines 



Physical properties from molecular 



This work makes heavy use of standard methods to derive phys- 
ical properties from molecular line data. Here we summarise the 
methodology used and gather the equations and constants in one 
place. All calculations were performed in SI units. 

A. 1. Kinetic temperature from A/H 3 

Firstly, we calculate the optical depth of the NH3 (1,1) transition 
via a comparison of brightness temperature T B in the main and 
satellite lines: 



1 



T B ,s 



1 -e" 



(A.l) 



The subscripts 'm' and V indicate the main and satellite groups, 
respectively, and 'a' is the theoretical T B S /T B m ratio under op- 
tically thin conditions. For NH3(1,1), the fractional intensity in 
the main group is 0.502 and in individual inner and outer satel- 
lite lines is 0.139 and 0.111, respectively. In practise we used 
the hyperfine fitting utility of the CLASS software package to 
simultaneously fit the five groups. It is often difficult to measure 
the optical depth directly from the NH3 (2,2) line profile as the 
signal-to-noise on the satellites is usually poor. Instead, T2,2 may 
be calculated from the brightness temperatures of the (2,2) and 
(1,1) main groups: 



T 2,2 



12,2 



I B 2 .2 main 
r Bl ! main 



(l-e-^u) 



(A.2) 



where f 1,1 = 0.502 and f2,2 = 0.796 are the fractional intensities in 
the (1,1) and (2,2) main groups, respectively, compared to their 
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satellites. The rotational temperature is found from the optical 
depth ratio T2,2It\,\ via: 



T|111 = -41K/lii||-^^ 



(A.3) 



where tj,k and Avj k are the total optical depth and the full-width 
half-maximum (FWHM) linewidth of the NH3 (J,K) transition. 
Finally, the rotational temperature may be converted directly to 
a kinetic temperature T^n via: 



I kin 



[1 + (T k i„/41 K)ln(l + C23/C21)] ■ {AA) 

Dan bv et"aTI (Q988) have used large velocity gradient (LVG) 
models to calculate the collisional coefficients C23 and C21, 
for kinetic temperatures ranging from 5-300K. To derive 
Equation IA.4I we assume that only the J,K = (1,1), (2,1) and 
(2,2) transitions are involved. At low temperatures this as- 
sumption is valid, however, above T lot ~ 30 K the excitation 
of NH3 to higher energy levels becomes non-negligible. This 
means that the above T rot ~ 30 K the kinetic temperature found 
via c omparison of the NH 3 (1,1) and (2,2) lines becomes unreli- 
able. Tafa lla et~aT] (|2004) have investigated the accuracy of this 
method and find that kinetic temperatures derived are accurate 
to better than 5 percent below 20 K. 

A.2. Column density and core mass 

If the optical depth and excitation temperature in any (J,K) tran- 
sition are known, the total column density of emitting molecules 
may be calculated assuming local thermal equilibrium (LTE) via 
the equation 



N 



8k7TV 



J,K 



A„ihc 3 



|k r Thdv ( I \ J__ 

c 3 J-oo b U -e- T / gj, K 



(E,, K )/kT 



gj,K Q(T rot ) 



(A.5) 



where Q(T rot ) is the partition function. The exact partition func- 
tion is the sum over all the levels however, this can be ap- 
proximated by an e xponential fit to the discrete values in the 



proximatea p y an e xponential tit to tne discrete values in 1 
iPickett. etaf] (U998) spectral line catalogue. For NH 3 this is: 



Q(T) = 0.1266 T 



A.3. COandHCO + 



(A.6) 



Physical conditions were calculated on a pixel-by-pixel basis in 
the Mopra data. Using the regrid task in miriad, the CO data- 
cubes were first resampled to the same pixel-scale in RA, Dec 
and velocity. The maths task was used to perform all subse- 
quent calculations. For each spatial pixel we calculated the opti- 
cal depth as a funct ion of frequency using a modified version of 
the isotope ratio dBourkeet aill 19971) : 

T MB (C 18 0) l- e - (T .* /x) 



T MB ( 13 CO) 



1 



(A.7) 



In the above equation, Tmb is the peak main beam brightness 
temperature and X is the abundance ratio of [ 13 CO] / [C 18 0] as- 
sumed to be constant. During the calculation pixels with insuffi- 
cient emission were blanked, resulting in a map of 13 CO optical 
depth with real values wherever C I8 emission was detected. 
The expect ed [ 13 CO/C 18 Q] ratio in giant molecular clouds is 
~ 1 1 .5 (see Goldsmith e tldll 19971 and references therein), how- 
ever, we found that this value resulted in negative optical depths 
in some parts of the map. The lowest abundance ratio that did 



not result in negative optical depths was ~ 16. This value was 
used in our calculations and should be considered a lower limit. 
Optical depths ranged from 0.1 < t < 4.4. Using the 13 CO op- 
tical depth and peak brightness temperature maps as input, we 
then constructed an excitation temperature map via: 



hv u 



In 1 



(hvu/k) 



T 13 /C T +J v (T bg ) 



(A.8) 



where C T = (1 - e~ T "). At this stage we blanked pixels with spu- 
riously excitation temperatures (T ex <4 K and T ex >40 K) due 
to noisy C 18 data which had passed our initial signal-to-noise 
cutoff. Assuming a thermalised population with T^n = T ex , a 
column density map was be produced via Equation IA.5I The 
partition function Q( T) for 13 CO has bee n measured exactly 
at discrete values by [Pickett, et alj dl998l) spectral line cata- 
logue. Intermediate values were determined via an exponential 
fit, which for CO is: 



Q(T kin ) = 0.3968 T, 



kin- 



(A.9) 



The CO column density map is presented in Figure [TCfl-a. 

We calculated the column-density of HCO + using the same 
procedure, and a partition function of the form: 



Q(T kin ) = 0.4798 

A kin- 



(A. 10) 



The [HCO + ] / [H 13 CO + ] was set to 40 for these calculations. 



A.4. N 2 H + and CS 

A slightly different approach was taken for the N2H + data. In a 
similar manner to NH3, the hyperfine structure in the spectrum 
of N2H + (1-0) allows the direct determination of optical depth. 
We fit all spectra in the N2H + data-cube using the HFS method 
in class and found that the N2FU emission was consistent with 
it being optically thin (r <k 0.1) over the entire cloud. Hence, 
no independent estimate could be made of the excitation tem- 
perature. N2FU and NH3 have similar critical densities and, be- 
ing nitrogen bearing molecules, are often detected under similar 
conditions. We calculated the column density of N2I-U assum- 
ing excitation temperatures equal to the NH3 rotational temper- 
atures. In practise, we smoothed the T rot (NH3) map to the same 
resolution as the Mopra N2H + data. We then produced a col- 
umn density map by combining this T ex map with the N2H + in- 
tegrated map via Equation [A3] in the limit that r <K 0. The par- 
tition function for N2H + w as again interpolated from the values 
in the lPickett. et alJd 19981) catalogue: 



Q(T kin ) = 4.198 T 



kin- 



(A.ll) 



A column density map of CS was prepared in a similar man- 
ner, although we note that our assumption of optically thin emis- 
sion may not be valid. Judging by the non-Gaussian shapes of 
some line profiles, regions with strong emission may have sig- 
nificant optical depths. 
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Table A.l. Molecular constants a . 



Species 


Transition 


Frequency 


EJk 


Aui 


n C rit'' 




(Jk) 


(GHz) 


(K) 


(s" 1 ) 


(cm~ 3 ) 


HCO + 


1 — > 


89 1 88526 


4.280 


3 020x 1 0~ 5 


3.0x10 s f 


ri CU 


1 — > U 


OA HZ. a iir\ 

oo. / j43 31) 


a 1 /in 
4. 16U 


Z.oUUX 1U 


3.UX1U ' 


13 CO 


1 -> 


110.210353 


5.288 


6.389xl0" 8 


6.5 X 10 2 + 


C 15 


1 — > 


109.782173 


5.267 


4.961x10 




cs 


2 -> 1 


97.980950 


7.052 


1.703xl0~ 5 


2.0 x 10 5 


N 2 H + 


li ->o, 


93.171880 


4.471 


3.654xl0~ 5 


2xl0 5 "" 




1 2 -> 0] 


93.173700 










lo-»Oi 


93.176130 








NH 3 


ii -» 1;, 


23.694480 


17.215 


1.67xl0~ 7 


2xl0 4 ™ 




2 2 -> 2',, 


23.722630 


63.866 


2.29xl0~ 7 





" Unless otherwise noted, all values quoted here are sourced fr om the NASA Jet 

Propulsion Lab (JPL) spectral line database (Pickett, et lil]|l998l) . 

^ Critical density values marked with a f are calculated from the formula 

Dcrit = A ul /y Hi , where y H is the 'standard' collisional rate coefficient for H 2 , assumed 

e qual to 10 10 cm 3 .s~ 1 .H2-m ol ecules~' . Oth er values are referenced as follows: 

"" lUngerechts et al.Ul997h . " 1Swadelll989l . 



